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Abstract:  From  an  analysis  of  the  effect  of  a  rough¬ 
surfaced  target  and  the  effect  of  atmospheric  turbulence 
on  a  propagating  laser  beam,  and  the  double-stochastic 
field  model,  we  propose  an  intensity  probability  density 
function  and  the  photoelectronic  count  statistical 
relationship  of  a  laser  beam  reflected  from  a  location 
target  and  propagating  through  a  weak  fluctuating 
turbulent  atmosphere . 

1 .  Introduction 


When  a  laser  beam,  reflected  from  a  location  target, 
propagates  through  a  turbulent  atmosphere,  its  power,  phase  and 
polarization  begin  varying  randomly,  due  to  stochastic 
fluctuations  in  the  atmospheric  refractivity  ratio  and  stochastic 
fluctuations  of  scattering  surfaces.  As  a  result,  this  process 
leads  to  several  effects  such  as  beam  expansion,  flickering, 
reduced  timewise  and  spatial  coherence,  and  whole-beam 
vibrations;  these  affect  the  capabilities  of  laser  radar,  namely: 
accurate  target  positioning,  image  resolution,  and  so  on,  and 


*This  paper  was  presented  at  The  2nd  National  Laser  Atmospheric 
Propagation  and  Engineering  Application  Seminar. 
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thus  have  attracted  growing  interest  [1-3].  We  discussed  the 
phase-fluctuation  statistical  distribution  [1]  of  backscattered 
waves  reflected  from  a  surface,  given  the  condition  that  the 
rough  surface  is  considered  as  a  .simple  approximation  of 
turbulent  phase  screening;  additionally,  we  analyzed  the 
statistical  distribution  law  [4]  of  the  backs cattered-wave 
optical  field,  given  the  condition  that  the  spatial  nonuniformity 
of  backscattered  wave  is  taken  into  account.  When  we  take  into 
account,  respectively,  laser  beam  modulation  by  atmospheric 
turbulence  and  by  the  location  target,  this  paper  gives  a 
theoretical  analysis  of  the  intensity-fluctuation  statistical 
distribution  law  of  backscattered  waves  reflected  by  a  location 
target  as  a  laser  beam  propagates  through  the  turbulent 
atmosphere,  allowing  for  laser  beam  modulation  by  atmospheric 
turbulence  and  by  location  target,  respectively. 

2.  Intensity  Probability  Density  Function  of  Circular  Gaussian 
Scattering  Spots 

When  a  laser  beam  propagating  through  a  turbulent  atmosphere 
is  reflected  by  a  location  target  surface,  it  will  be  subjected 
to  interference  owing  to  atmospheric  turbulence  and  surface 
fluctuations.  Since  the  target-surface  fluctuation  law  and  the 
atmospheric-turbulence  stochastic  fluctuation  law  are 
statistically  independent,  we  calculate  their  interference  on  the 
laser  beam  separately.  During  the  analysis,  we  regard 
atmospheric  turbulence  prior  to  and  after  scattering  by  the 
target  as  the  same  interference  factor  and  take  the  effect  of 
target  stochastic  fluctuation  surface  on  scattering  waves,  as 
another. 

Below  we  give  an  analysis  of  the  statistical  law  of  laser- 
beam  backscattered  waves  from  a  location  target  surface  that  is 
made  up  of  an  optically  rough  part  and  an  optically  smooth  part. 

As  an  approximation,  the  analysis  is  limited  to  surface 
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scattering  in  which  the  backscattered  wave  spots  satisfy  the 
circular  gaussian  distribution  without  considering  shadowing 
effects  due  to  the  very  long  distance  between  observation  surface 
and  scattering  surface.  The  backscattered  wave  total 
instantaneous  field  E-arrow  is  a  complex  quantity  consisting  of 
backscattered  waves  reflected  by  the  rough-part  of  the  target 
surface  and  by  the  smooth  part  of  the  target  surface,  namely, 

E  —  +  E,  —  £exp(j0)  ^  i  ^ 

where  E  and  B  are,  respectively,  the  amplitude  and  the  phase  of 
E,  Ed  and  Eg,  respectively,  are  the  scattering  stochastic 
component  of  the  rough  part  and  the  backs cattered-wave  radial 
component  from  the  smooth  part. 

Let  the  stochastic  scattered  component  Ed  in  E  be  a  compound 
gaussian  stochastic  quantity  of  zero  mean  value  within  the 
observation  surface;  the  real  Erd  and  the  imaginary  E^  part  of 
Ed  are  gaussian  variables,  mutually  independent  and  evenly 
distributed  in  space  while  the  radial  component  jf  consists  of  N 
backscattered  waves  from  the  optically  smooth  part,  then 

E>  =  2  ^.exP0»)  =  Etexp(jtp)  ^  2 

where  An  and  Phin  are  the  fixed  constant-amplitude  and  phase  of 
the  n-th  reflected  laser  beam  from  the  smooth  part,  while  Eg,  phi 
are  the  amplitude  and  phase  of  the  spatial  stochastic 
distribution. 


According  to  [5],  the  light  intensity  of  field  i  obeys  the 
Rice  distribution 


p*cn 


</<> 


exp  ( 


-  i+JiWh/K 


</,> 


<f<> 


(3) 


where  /  =  El,lt  =  £,2, </,)  =  are  the  variance  of  stochastic 

scattering  component  Ed,  l0(  »)  is  the  Bessel  function  of  first 
category  zero-order  deformation. 
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Apart  from  the  foregoing  description  of  beam  scattering  by  a 
location  target,  the  effect  of  atmospheric  turbulence  on  a 
propagating  laser  beam  should  also  be  calculated.  To  simplify 
analysis,  let  us  assume:  1)  the  beam  propagates  only  in  a  weakly 
turbulent  atmosphere,  during  which  modulation  of  the  time- 
averaged  intensity  value  by  beam  vibration  can  be  ignored;  2)  the 
linearity  of  beam  emitted  to  the  target  is  smaller  than  target 
linearity;  and  3)  the  spatial  distribution  of  target  fluctuation 
is  even.  According  to  the  foregoing  approximation,  Ed  is  a 
gaussian  variable  of  zero  mean  value,  and  it  is  assumed  to  be 
known  that  the  effect  of  atmospheric  turbulence  on  backscattered 
wave  from  the  target  surface  is  virtually  an  effect  on  the 
reflection  component  in  the  reflector  direction  Eg,  in  other 
words,  Eg  is  assumed  to  be  a  time-varying  stochastic  variable 
modulated  by  atmospheric  turbulence,  and  Eq.  (3)  is  the 
conditional  Rice  distribution,  then 


?o  (///,) 


1 

</<> 


•exp(  — 


l±L 


)/«( 


2V/T 

</<>  ; 


(4) 


While  the  unconditional  probability  density  function  of  the 
backward  wave  intensity  is 

?(/)  =  (5> 

where  pt(Ig)  is  the  Ig  fluctuation  probability  density  function 
under  interference  by  atmospheric  turbulence. 


On  the  assumption  of  a  split  scattering  model,  the  light 
wave  radiation  field  at  a  given  point  and  at  a  given  instant 
within  the  observation  surface  can  be  considered  as  being  formed 
through  the  superimposition  of  scattering  along  different 
turbulent  paths,  while  the  fluctuation  distribution  of  total 
scattering  field  intensity  x  can  be  expressed  with  the  deformed 
Rice-Nakagami  distribution  [6]  as  follows: 
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1  r  C4*  +  r)-|  ,  ,2AJT' 

p,(x)  =  —exp  [• - - - -J/oC— —  )  (6) 

where  A  is  nonscattered  coherent  constant-amplitude  component  and 
b  is  the  mean  wave  intensity  of  the  circular  compound  gaussian 
stochastic  component  in  the  relation  as  <*>  =  A2  +  b,a,2  =  a2  4-  2Azb, 


According  to  N.  S.  R.  Gudimetla  et  al.  [7],  the  14- 
distribution  is  a  desired  approximation  for  the  deformed 
Rice-Nakagami  distribution.  Thus,  the  M-distribution  can  be  used 
to  approximate  the  Rice-Nakagami  distribution  as  follows : 

.  _  Af*xv-,ex p(—  Mx  /  (x  >) 

r,M - FTr)<«>” —  (?) 


In  this  way,  the  fluctuation  distribution  of  target  radial 
reflection  component  intensity  Ig  can  be  derived  as  follows: 


pXI  )  =  g**r*~'«xp(-  Ml./U.)) 

’  rm  )</,>*. 


(8) 


The  probability  density  function  of  the  backscattered  wave 
intensity  thus  derived  is: 


2  •  II 
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(9) 

where  A/t.„(z)  is  the  Whittake  function  and  r(z)  is  the  gamma 
function. 

The  nth  factorial  moment  of  the  backward  wave  intensity  is: 


5 


</*>  — 
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where  2Fi(a/  b,  c,  x)  is  the  hypergeoxnetric  function. 
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If  the  location  target  is  much  smaller  than  transmitting 
beam,  then,  taking  advantage  of  the  statistical  independence  of 
the  atmospheric  refractivity  ratio  fluctuations  and  the  target 
surface  fluctuations,  light  wave  scattering  by  target  can  be 
taken  as  re-modulation  of  the  radial  component  of  scattered  wave 
in  atmospheric  turbulence.  Similarly,  by  writing.  Eq.  (6)  as 
Pt=Pi(x/A2)  and  by  writing  Eq.  (4)  as  Po(A2)  in  the  M-distribution 
approximation,  the  intensity  probability  density  function  of 
backs cattered-wave  field  can  be  obtained  as  follows: 


/  1  \*  /  7  \  T  -  7 <7,  >  1 

■  p2<7>  =  ^"(<7<>  +  AfV/,))  <z lexH-w)  ,f,r  >(</<> +  m<7,>)J " 

where  ^(a^^)  is  the  confluent  hypergeometric  function. 


(11) 


Individual  factorial  moments  of  light  intensity 
corresponding  to  Eq.  (11)  are: 

</•>  -  {a  +  VhF'  ["•’ +  +  «(/,>]  <12> 

Since  our  analysis  focuses  on  the  intensity  fluctuation  of 
backscattered  wave  from  a  location  target  during  propagation 
through  a  stable  atmosphere  with  turbulence  from  gentle  to  weak, 
the  integral  intensity  corresponding  to  Eq.  (9)  has  the  same 
distribution  as  instantaneous  intensity  [8],  W=  J +  7(U  ,  that 


P.GD 


According  to  Mandel's  equation,  the  optoelectronic  count 
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distribution  is 


Pi(ot,T)=["  — e~'*rvi(W',<lW 

Jo  ml  '  1 


Mmd 


P  (  TTl  — p  1  )  i 


m9  "1“  Mnif  J  m  ]  mg  ^  1  +  m* 


i  +  1 


:ZF,  m , 


m  +  1,1,-= 


(m,  4*  Afm^Cl  +  m 


where  a  is  the  ratio  factor,  m*-a<-w*'>>m'-a<w,>' 

In  the  same  way,  p2(m,T),  which  corresponds  to  Eq.  (11),  can  be 
derived . 


3.  Intensity  Probability  Density  Function  of  Elliptic  Gaussian 
Scattering  Spots 


Generally  speaking,  the  compound  amplitude  of  scattered 
spots  from  a  rough  surface  can  satisfy  the  elliptic  gaussian 
distribution  [9]  and  therefore,  it  is  necessary  to  analyze  the 
intensity  probability  density  function  of  backscattered  waves 
reflected  from  by  this  target  during  beam  propagation  through  a 
turbulent  atmosphere.  Since  we  can  somehow  simplify,  through 
coordinate  conversion,  the  general  elliptic  gaussian  distribution 
into  the  product  of  distribution  functions  corresponding  to 
independent  components,  the  combined  probability  density  function 
of  the  real  part  and  the  imaginary  part  of  scattered  spot 
amplitude  can  be  expressed  as: 


2^"exP  [- 


(AF,)* 


2<r,2  J 


p(E,Jp(E,J 


(14) 


where  Ez  and  Ey,  respectively,  are  the  real  part  and  the 
imaginary  part  of  E,  =  E,  —  < Et >  is  the  fluctuation 

quantity  of  Ez,  <f,>  =  .Vexp(— a2/2)  is  mean  value  of  the  Ez 
component,  =  (IV / 2)[1  4-  exp(—  2a2)  —  2exp(—  (j2)]  is  fluctuation 
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variance  of  Ez,  a-,2  =  <^,2>  =  (A/2)[l  —  exp(— 2c^)]  is  fluctuation 

variance  of  Ey,  <*>  is  the  system  statistical  average,  N  is  the 
number  of  target  scattering  elements,  o  is  the  phase 
fluctuation  standard  difference  of  target  scattering  field.  By 
making  an  analysis  similar  to  that  described  in  reference  [4],  we 
acquire:  . 


,(1,0)  =^-1— /exp  L  2(/  +  7,-2  7^cos(g-^)) 
l  (1+ ?>)</,> 


—  exp 


_  2(7  +  /,  -  2  cos(fl  -  ftyy  \ 
(1-7>X7<>  JJ 


(15) 


By  using  relations  between  joint  probability  density  and  marginal 
probability  density,  the  monovalent  probability  density  function 
of  light  intensity  can  be  derived  as  follows: 


p(/) 


.(1  +9>)<7<>J 


1  /  r  2(7  +  7.)  1  , 

V  <^  >  \exp  L  (1  +  9>)<7<>J/|1 

_  exp  r_  2(7  +  V.1  /  4  VtTt  1  \ 

L  ( 1  -  9>)  <7,  >J  [( I  _  ,,)</<  )J  / 


(16) 


(1  -<P)<7, 

Upon  consideration  of  Ig  modulation  by  atmospheric  turbulence, 
when  the  target  is  larger  than  beam  linearity,  the  intensity 
probability  density  function  of  backscattered  waves  by  a  location 
target  in  turbulent  atmosphere  propagation  can  be  derived  from 
Eq.  (5)  as  follows: 
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Correspondingly,  the  various  intensity  factorial  moments  are 
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(18) 


Similar  methods  can  also  be  used  to  obtain  corresponding 
results  when  the  location  target  is  smaller  than  beam  linearity. 


4.  Conclusion  and  Discussion 


Based  on  the  statistical  independence  of  atmospheric 
turbulence  fluctuations  and  target  surface  fluctuations,  as  well 
as  their  separate  modulation  capability  over  a  propagating  laser 
beam,  this  paper  has  derived  the  flicker  probability  density 
function  of  backscattered  waves  reflected  by  a  rough  surface 
during  propagation  through  turbulent  atmosphere  as  expressed  in 
Eqs .  (9),  (11)  and  (17).  These  results  indicate  that  intensity 
fluctuation  distribution  under  joint  interference  by  atmospheric 
turbulence  and  location  target  proves  to  differ  from  the 
propagating  laser  intensity  fluctuation  distribution  under  the 
single-existence  of  either  stochastic  rough  surface  or 
atmospheric  turbulence.  This  deserves  attention  when  strict 
calculations  are  conducted  over  the  intensity  density  function  of 
backscattered  waves  reflected  from  a  target  during  propagation 
through  turbulent  atmosphere . 

Our  analysis  shows  that  beam  vibration  affects  the 
measurement  of  average  light  intensity  and  therefore  affects  the 
light  flicker  probability  distribution  [6];  also,  using  the  M- 
distribution  as  replacement  of  the  approximate  Rice-Nakagami 
distribution,  given  the  condition  that  the  spherical  wave  source 
and  parameter  M>5  appears  to  be  a  desired  approximation  with  an 
error<3%.  However,  in  our  analysis,  the  beam  vibration  effect  is 
ignored,  and  the  M-distribution  is  adopted  to  approximate 
Rice-Nakagami  distribution,  the  results  thus  obtained  are 
applicable  only  to  flickering  probability  distribution  of  the 
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reflected  beam  from  a  rough  surface  during  divergent  spherical 
wave  propagation  through  a  weakly  turbulent  fluctuation  region. 
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